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基于 Givens 变换 和 二 阶 振荡 W-C-PSO 优化 的 盲 源 分 离 算 法 
张 华 伟 ， 张 天 器 ， 间 董 华 


(重庆 邮电 大 学 信号 与 信息 处 理 重庆 市 重点 实验 室 , 重庆 400065) 


摘 要 : 针对 智能 算法 在 实现 讶 源 分 离 时 容易 陷入 局 部 最 优 ， 且 收 伊 速度 缓慢 的 问题 ， 提 出 了 一 种 基于 Givens 变换 和 
二 阶 振荡 粒子 群 优化 的 育 源 分 离 算 ie sw sw 
迁 代 ， 来 提高 算法 的 整体 性 与 全 局 搜索 能 力 ; 再 引入 二 阶 振荡 环节 增加 种 群 的 多 样 性 ， 这 样 算法 不 易 陷 入 局 部 最 优 ; 
此 外 ,采用 Givens 变换 将 分 离 矩 阵 转换 成 旋转 角度 表示 形式 ， 来 降低 算法 的 复杂 度 。 仿 真 表明 ， 该 算法 能 有 效 实现 机 
械 振动 信号 和 语音 信号 的 讶 分 离 ， 并 且 相 比较 其 他 算法 具有 更 快 的 收 化 速度 和 更 好 的 分 离 性 能 。 

关键 词 : 讶 源 分 离 ; 粒子 群 算法 ; 二 阶 振荡 ; Givens 变换 

中 图 分 类 号 : TN911.7 doi: 10.19734/i.issn.1001-3695.2018.05.0447 


Blind separation method based on based on Givens transformation and second-order oscillatory 
W-C-PSO 


a Zhang Huawei, Zhang Tianqi, Liu Donghua 
(Chongqing Key Laboratory of Signal & Information Processing, Chongqing University of Posts& Telecommunications, 
Chongqing 400065, China) 


Abstract: For the intelligent algorithm, it is easy to fall into local optimum when implementing blind source separation, and 
the convergence Speed is slow, this paper proposes a blind source separation algorithm based on Givens transform and 


second-order oscillator particle swarm optimization. This algorithm constructs the functional relationship between the inertia 


i weight and the learning factor, and adjusts the algorithm together to improve the algorithm's overall and global search ability. 
The second-order oscillation increases the diversity of the population, so the algorithm is not easy to fall into the local 


optimum. In addition, the algorithm uses Givens transformation to convert the separation matrix into a rotation angle 


representation to reduce the complexity and accelerate the convergence speed. Simulation results show that the algorithm can 
effectively achieve the blind separation of mechanical vibration signals and speech signals, and has faster convergence speed 
and better separation performance than other algorithms 


Key words: blind source separation; particle swarm optimization; second-order oscillatory particle swarm; givens 


transformation 
0 引言 虑 到 根据 信 源 统计 特性 选取 非 线性 函数 的 问题 ,影响 因素 较 多 ， 
加 并 与 源 信号 信息 未 知 相 矛 盾 。 


盲 源 分 离 是 指 在 源 信号 和 传输 信道 未 知 的 条 件 下 ， 从 混合 当前 更 多 的 研究 者 将 智能 算法 应 用 到 盲 源 分 离 中 来 解决 上 
言 号 中 分 离 出 各 路 待 估计 信 源 的 过 程 ， 在 参数 估计 、 各 类 信号 述 问题 651。 文献 [7] 提 出 了 基本 粒子 群 (particle swarm 
分 析 、 故 障 检测 等 领域 均 具 有 重要 应 用 二 9。 当前 传统 的 盲 源 optimization, PSO) 算法 ， 较 蚁 群 算法 、 蜂 群 算法 等 智能 算法 ， 
分 离 算法 研究 的 较为 成 熟 ， 如 快速 独立 分 量 分 析 法 、 自 然 梯度 PSO 算法 的 收敛 速度 更 快 ， 全 局 搜索 能 力 更 强 ， 和 迭代 过 程 更 简 
法 、 等 变 自 适 应 算法 等 。 此 类 算法 使 用 自 适 应 变 步 长 改进 了 算 ”单一 些 。 文 献 [9] 提 出 了 基于 自 适 应 调整 惯性 权重 的 PSO 算法 ， 
法 的 收敛 速度 申 ， 增 加 动量 项 提高 了 算法 的 稳 态 误差 由 ,但 考 。 虽然 解决 了 非 线性 函数 选取 问题 ， 但 该 算法 的 收敛 性 并 没有 很 
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大 改善 且 容 易 陷入 局 部 极 值 ; 文献 [9] 采 用 自 适 应 遗传 变异 环节 
增加 粒子 群 的 多 样 性 ， 提 高 了 算法 的 分 离 性 能 ， 但 是 
算 变 异 概率 和 子 代 的 关系 ， 算 法 的 复杂 度 也 急速 上 升 。 
针对 上 述 算法 的 缺点 , 本 文 提出 了 一 种 基于 Givens 变换 和 
二 阶 振荡 W-C-PSO 优化 的 盲 源 分 离 算 法 。 该 方法 采用 粒子 群 
算法 对 目标 函数 进行 优化 ， 将 惯性 权重 与 学 习 因 子 非 线性 函数 
结合 ， 共 同调 节 粒 子 的 迭代 ， 这 样 既 可 以 提高 算法 的 整体 性 ， 
又 可 以 增强 全 局 收敛 能 力 ; 其次， 引入 二 阶 振荡 环节 增强 粒子 
的 多 样 性 , 使 粒子 极 易 跳出 局 部 最 优 ; 并 采用 Givens 变换 将 分 
离 矩阵 使 用 旋转 角度 表示 ， 来 降低 算法 的 复杂 度 ， 提 高 算法 的 
收敛 速度 。 最 后 ， 采 用 本 文 算法 对 振动 信号 和 语音 信号 进行 盲 
分 离 ， 实 验证 明 该 算法 较 其 他 算法 具有 较 好 的 收敛 速度 和 全 局 
搜索 能 力 ， 并 有 效 解决 易 陷入 局 部 极 值 的 问题 


1 ” 盲 源 分 离 模 型 


orn 省 述 为 : 假设 有 个 信号 源 和 mm 个 传感器 ， 则 
观察 信号 与 源 信号 呈现 的 表达 式 为 


x(1)=As(1) +n(?) (1) 


其 中 : s(7?)=[si,s,,…， 


s,] 为 mxl 维 源 信号 向 司 


x(1) =[%,%,..., 


和 矩阵 ，n(z) 为 背景 噪声 ， 在 这 里 本 文 不 予 考虑 。 在 信道 性 质 和 
源 信号 均 不 知 时 ， 仅 由 观测 数据 x(7) 和 相应 的 学 习 算 法 得 到 最 


元 ] 为 mxl 维 观测 数据 向 量 ; A 为 mxn 维 混合 


佳 分 离 矩阵 W ， 使 得 算法 输出 的 (人) 是 源 信号 的 最 佳 估计 。 
y(t)= Wx(1) Q) 
其 中 ，y(?)=[y,y,…,y] 。 所 以 ， 盲 源 分 离 过 程 转换 成 了 在 


某 一 个 独立 性 判 据 意义 下 , 通过 相应 的 算法 对 分 离 矩 阵 W 进行 
寻 优 的 过 程 ， 使 各 分 量 实现 相互 独立 。 本 文采 用 负 业 9 作为 分 
离 出 的 估计 信号 相互 独立 的 依据 来 建立 目标 函数 。 


一. 


Ey)= ek (Cy) G) 
2 ) 
其 中 的 霄 度 总 (7 到 3 在 E(yy )=1 条 件 下 ， 
yi 


目标 函数 f(y) 越 大 ， 则 各 分 离 信号 y 之 间 的 统计 独立 性 越 高 ， 
也 就 是 分 离 结果 越 好 。 


2 基于 二 阶 振荡 改进 的 W-C-PSO 算法 
2.1 基本 粒子 群 算法 


PSO 算法 根据 目标 函数 计算 适应 度 值 ， 使 粒子 不 断 地 更 新 
速度 和 位 置 ， 最 终 找到 全 局 最 优 位 置 。 其 速度 和 位 置 的 更 新 公 
式 如 下 : 
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V(t+1) =V,,(D)+on|p;, —X,,(n)|+ 
Cn [Lp 本 X,,(7)| 

X(t+1)=X, (t+V(t+1),j=1,2,.…,D (5) 
其 中 : oc 和 cc 分 别 为 认 知 学 习 因 子 和 社会 学 习 因 子 ，7 和 是 

分 布 于 [0,1] 的 随机 数 ; D 为 粒子 搜索 空间 维度 
X= (Xii, Xs,…,Xip ) 为 第 i 个 粒子 位 置 ， 郧 =(VV Von) 
为 第 i 个 粒子 速度 ，P =( pi pio,…, Pip) 为 粒子 个 体 最 优 位 置 
P, =(poi,Pes…, Pen) 为 当前 种 群 全 局 最 优 位 置 。 
传统 PSO 算法 中 令 学 习 因 子 ac 和 c 为 常数 ; 文献 [11] 将 1 


性 权重 @ 应 用 到 PSO 算法 中 ,有效 平衡 了 算法 的 收敛 速度 和 
索 能 力 。 


(4) 


六 洒 和 


V(t +1)= @() V(t) + on |p;, —X,) 人] 
+on| ps, -xX,,(7)| 
其 中 : 中 采用 下 式 线性 递减 的 形式 ， 即 在 [0.95,0.4| 内 变化 。@ 
大 变 小 ， 使 算法 在 前 期 扩大 搜索 范围 ， 全 局 搜索 能 力 较 好 ; 
而 在 后 期 集中 局 部 搜索 ， 提 高 准确 性 。 


(6) 


> 


1 
ow( 1) Cnax ( Cnax hin ) X N (7) 


其 中 ， + 为 迭代 次 数 ，WN 为 最 大 友 代 次 数 : wu 为 惯性 
大 值 ， wmn 惯性 权重 最 小 值 。 
2.2 改进 粒子 群 算法 

可 知 惯性 权重 可 有 效 改 善 搜索 能 力 , 而 由 式 (4) 可 知 学 习 
因子 6 、2 分 别 为 调节 粒子 向 自身 最 优 方 向 、 全 局 最 优 位置 的 
前 进步 长 ,因此 调节 学 习 因 子 可 以 有 效 避 免 粒 子 陷 入 局 部 极 值 ， 
改善 算法 搜索 性 能 。 但 在 以 往 的 算法 改进 中 ， 这 两 个 参数 互 不 
关联 ， 这 在 某 种 角度 上 说 降低 了 PSO 算法 的 整体 性 ， 不 利于 粒 
子 的 全 局 搜索 ， 因 此 本 文 提出 将 这 两 个 参数 非 线 性 结合 来 增强 
算法 的 整体 性 能 

| al(I)=Yo (7) +no(t) +e 


Cs (7) =2.5 一 Cl (@) 


区 
于 
沽 


(8) 


其 中 : Y、n 、& 为 常 系数 。 通 过 测试 本 文 惯 性 权重 采用 下 式 
指数 递减 形式 : 
(1) = On, + (rar -omj*exp| -20(117) | (9) 


这 就 是 本 文 所 说 的 W-C-PSO 算法 。 

在 粒子 群 算法 中 ， 随 着 粒子 的 不 断 搜索 ， 粒 子 的 种 群 多 样 
性 在 不 断 减 小 ， 这 将 不 利于 继续 搜索 。 因 此 本 文 在 W-C-PSO 
算法 的 基础 上 引入 二 阶 振荡 环节 ， 可 以 使 粒子 在 下 一 次 迭代 中 
向 着 较 优 的 位 置 移动 ， 即 增强 粒子 的 多 样 性 ， 当 搜索 范围 扩大 
时 , 二 阶 振荡 有 助 于 粒子 跳出 局 部 , 改善 算法 的 全 局 搜索 能 力 。 
改进 后 粒子 群 速度 迭代 可 表示 为 


V(r+1)=0()V()+ 
)r [pi (+E)X, (+X GD (10) 
(Da[ psy-(1+6)X,,() + EX,, (il)| 
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其 中 : 所 ， 允 为 随机 数 。 在 该 算法 欠 代 前 期 〈 小 于 最 大 友 代 次 求解 的 参数 个 数 减少 了 
数 的 2 ) 选取 号 < 人 5 ， 吉 < 纪 25-1 ， 可 提高 算法 的 二 二 (10) 
Cn Cp 2 2 
全 局 搜索 能 力 ; 在 迭代 后 期 (大 于 等 于 最 大 迭代 次 数 的 2 ) 这 样 不 仅 保证 了 W(?) 的 正 交 性 ， 而 且 使 其 求解 的 未 知 参 
2cr 一 1 Cr 数 减 半 ， 降 低 了 复杂 度 。 
选取 > 一 一 ， ee ne 22 ~， 可 加 快 算法 的 收敛 速度 。 
3.2 基于 二 阶 振荡 W-C-PSO 优化 的 讶 源 分 离 
基于 二 阶 振荡 W-C-PSO 算法 用 粒子 位 置 表示 正 交 分 离 矩 
3 Givens 变换 和 二 阶 振荡 W-C-PSO 优化 的 盲 
I 阵 ，Givens 变换 后 其 参数 表示 使 粒子 搜索 维度 D 从 n? 降 到 了 
n(n 一 1)/2, 提高 了 算法 效率 ; 惯性 权重 与 学 习 因子 两 个 参数 的 
3.1 Givens 变换 结合 提高 了 算法 全 局 搜索 能 力 ;， 引入 二 阶 振荡 环节 提高 了 种 群 
式 (3) 可 知 , 处 理 n 阶 和 矩阵 时 需求 解 n? 个 未 知 数 , 计算 量 多 样 性 ， 使 其 不 易 陷入 局 部 极 值 。 算 法 可 分 别 有 效 实现 对 动态 
较 大 。 为 此 ， 对 混合 矩阵 采取 QR 分 解 。 令 白化 矩阵 为 童 号 和 语音 信号 的 言 分 离 ， 有 具体 流程 如 图 1 所 示 。 
z(t)=Px(t)=PAs(t)=Us(7) (11) (FR ) 
利用 QR 分 解 ， 则 z( 四 =Us(7)=QRs(1)。 其 中 , P 为 白化 
| 观测 信号 去 均值 预 白化 处 理 | 


和 矩阵，Q 为 酉 矩阵 ，R 为 上 三 


E[z(1)z(?) |=E|(QRs()(QRs()) | 
=E|QRs(i)s() R'Q 
=QRE|[s(i)s(?) |R'Q 

源 信号 之 间 不 相关 可 知 ，E| s(n)s(?) |= 7， 得 到 R 为 本 


矩阵 ， 且 为 对 角 阵 。 上 式 可 
Qz() 


表示 为 : 
=Rs(?) 


分 离 信号 与 
信号 的 估计 ， 即 分 离 信号 。 


(14) 所 示 。 


1 0 0 
0 cosO sing 
G; = 2 
0 —sinQ cos0. 
[| 0 0 


则 分 离 矩阵 W(7) 可 表示 为 
WO=G-(OcC-(0…G(OD 
G (0 二 Qi， (NG, (人 ‘Gi, (7) 


G,(1)=6G,,(1)G;, (7)-…G, (7) 


GG, (7) G1 (7?) 
其 中 : 9 e[0,2x] 。Givens 变换 后 ， 每 个 Givens 入 


j 阵 。 经 过 白化 后 有 


加 


待 求 参数 ，G 有 站-1 个 ， 


个 ， 所 以 分 离 矩 阵 W(7) 中 的 求解 参数 有 : 


_ n(n— 
(+2+…+(-2)+(-D)xl= 一 7 


] Givens 变换 ， WI(1) 中 的 求解 


G, 有 n-2 个 ， 依 次 类 推 


(12) 


(13) 


源 信 号 之 间 成 比例 关系 ， 则 Rs(z) 可 看 作对 源 
则 Qz 可 看 作 分 离 矩阵 的 估计 ， 即 
分 离 矩阵 为 正 交 和 矩阵。 任何 一 个 行列 式 不 等 于 1 的 n 阶 正 交 算 
阵 均 可 以 转换 成 Givens 旋转 矩阵 的 乘积 形式 ， 记 为 G12， 


如 式 
0 
: (14) 
0 
1 
E 阵 中 有 1 个 
>» Gi 有 1 
y (15) 


参数 有 nn? 个 ， 所 以 


采用 Givens 和 矩阵 表示 初始 粒子 位 置 | 


计算 适应 度 值 


V 
| 记录 粒子 的 个 体 极 值 和 全 局 极 值 


站 


en ~ 
一 生理 找到 最 优 解 一、 了 
“一 或 迭代 结 


i 


于 
外 -一 


自 适 应 更 新 惯性 权重 


随 之 更 新 认 知 学 习 因 子 和 社会 学 习 因子 


引入 二 阶 振荡 环节 增加 种 群 


多 样 性 防 . 
陷入 局 部 最 优 


上 粒子 


图 1 基于 Givens 变换 和 二 阶 振荡 W-C-PSO 算法 流程 图 


v 
输出 分 离 矩阵 


计算 得 到 分 离 信 号 


人 结束 ) 


Fig.l Flow of W-C-PSO algorithm based on givens transform and 2- 


order oscillation 


a) 对 观测 到 的 动态 混合 信号 x(7) 中 心 化 和 白化 操作 。 


b) 初 始 化 学 习 因 子 a 、 
wumsx 和 最 小 值 ww ， 随 机 设置 


Cc 以 及 惯性 权重 中 ， 并 设置 最 大 值 
些 数量 的 分 离 矩阵 作为 各 粒子 


的 初始 位 置 ， 


采用 Givens 变换 后 的 参数 表示 。 


9 根据 式 (3) 求 解 适应 度 


dj 对比 每 个 粒子 的 适应 度 值 与 个 体 极 值 P; 和 全 局 极 值 了 , ， 


取 优 更 新 Pi、 


6e) 当 t>2 时， 根据 式 (8) (9) 分 
阶 振荡 ， 根 据 式 (10) (5) 分 别 更 新 粒子 的 速度 和 位 


f) 引 入 二 


吧 进 入 迭代 循环 ， 满 足 终止 条 件 


8 


) 别 调节 8 和 c、c 


| 


输出 分 离 结果 ; 否则 返 


步骤 c) 继 续 搜 索 
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4 ”仿真 实验 


为 验证 本 文 提 出 的 基于 Givens 变换 和 二 阶 振荡 W-C-PSO 
算法 的 有 效 性 ， 对 机 械 振动 信号 和 语音 信号 分 别 进行 盲 分 离 实 


验 。 本 文采 用 时 变 混 合 矩 阵 ， 表 示 为 
1.0 0.5 
A=Aotalrandrsize(ho 彰 ， 设 初始 -| 沁 | 


a =0.003 ，rand 必 icel4o] 为 与 40 同 维 的 随机 和 矩阵， 也 就 是 混 
合 信号 是 随机 动态 变化 的 。 


4.1 振动 信号 的 盲 源 分 离 实验 
本 实验 将 模拟 机 械 振动 信号 作为 讶 分离 实验 的 源 信号 。 依 
据 振动 信号 模型 03 构 造 的 两 路 源 信号 为 
s(t)=3(1.5+1)cos(250xt —2) 
s (1) = sin(150z7)[1 + sin (20z1) | 1 
本 文 算法 参数 设置 如 下 : 粒子 种 群 规模 为 20， 惯 性 权重 最 
大 值 wu =0.9 ， 惯 性 权重 最 小 值 wu =0.4 ， 学 习 因 子 和 迭代 参 
数 Y=05，m=1l1，s=0.3， 代 进化 次 数 为 300。 实 验 结果 如 
图 2~4 所 示 。 
5| | pp PI | 
| | Wi ] 
贡 oh | | 
:于 
0 200 
1 | | | | | 外 | | 
由 | hoa 省 0.4 
到 
| 本 本 硬 呈 而 交 间 二 
0 0.1 9 0.4 0 200 
图 2 源 信号 时 域 波形 和 频谱 
Fig.2 Time domain waveform and spectrum of source signal 
加 | i | | 
上 | | 中 山 | 人 
0 01 02， 03 04 0 wo 200 
Tm 1 
: 玫 | | | 
0 0 -让 03 04 0 “一 200 
ts YHz 
图 3 混合 信号 时 域 波 形 和 频谱 
Fig.3 Time domain waveform and spectrum of composite signal 
比较 图 2 源 信 号 和 图 4 分 离 信号 ， 从 时 域 中 可 以 看 出 两 个 
信号 的 波形 基本 一 致 ， 从 频 域 中 也 可 以 看 到 信号 在 同样 频率 点 
具有 最 大 峰值 ， 说 明 本 文 算 法 可 有 效 实现 对 机 械 振动 信号 的 言 
分 离 。 
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下 | 
划 | 
时 “ I 
1 从 中 
100 200 
fHz 
于 
| 
世 0 
畦 | 
1 由 
WD 
0 01 02 03 0.4 0 100 200 
t/s /Hz 
图 4 分离 信号 时 域 波 形 和 频谱 
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考 第 工期 


Fig.4 Time domain waveform and Spectrum of discrete signal 


的 效果 ， 将 对 两 路 语音 信号 


4.2 ”语音 信号 的 盲 源 分 离 实验 
为 验证 本 文 算法 在 实际 应 用 
进行 盲 源 分 离 , 其 参数 设置 如 4.1 节 。 实验 结果 如 图 5~8 所 示 。 
源 信 号 1 
妆 下 
20 0.5 1 1.5 2 
采样 点 数 10% 
源 信 号 2 
2 i 
- 细 的 
20 0.5 和 1.5 2 
采样 点 数 x 10 7 
图 5 源 信号 波形 
Fig.5 Waveform of source signal 
混合 信号 1 
0|- ee | 
20 0.5 15 2 
采样 点 数 10° 
混合 信号 2 
1 Ti | 中 
-PR a 喘 | 
采样 点 数 


图 6 混合 信号 


Fig.6 Waveform of composite signal 


Fig.7 Mixed signal 


蜡 混 信号 1 ， | 

- 1 hi ~ | 

WI | | 由 

下 由 | l | | 

0 0.5 采样 点 数 15 
图 7 解 混 信号 
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图 8 分 组 散 点 图 
Fig.8 Block scatter plot 
比较 图 5 和 7， 从 波形 中 可 以 看 出 信 


号 成 功 分 离 出 来 。 并 


且 依 据 真实 听 到 的 音效 ， 幅 度 的 模糊 可 以 忽略 。 图 8 为 分 离 信 
号 的 分 组 散 点 图 ， 图 中 两 条 直线 说 明了 源 信号 和 分 离 信号 有 很 
高 的 相似 性 。 
4.3 算法 对 比 

1) 收 敛 性 与 搜索 能 

采用 PI 曲线 作为 性 能 指标 031， 其 定义 为 : 


jc 


-2 ] bc | 


其 中 :K 为 欠 代 次 数 ，G( 让 ) 为 k 次 迭代 时 全 局 矩阵 的 第 (i,j 
个 元 素 。 本 实验 将 对 比 传统 基本 PSO 算法 、W-C-PSO 算法 以 
及 引入 二 阶 振荡 的 本 文 算法 的 收敛 性 及 全 局 搜索 能 力 ， 如 图 9 
所 示 。 


(18) 


rt 


0 50 100 150 200 250 300 


迭代 次 数 
图 9 各 算法 收敛 曲线 1 
Fig.9 Convergence curves of algorithms 

从 图 中 可 以 看 到 , 基本 PSO 算法 的 曲线 没有 达到 最 优 值 就 
收 化 了 ， 导 臻 算法 陷入 局 部 极 值 ， 采 用 惯性 权重 与 学 习 因 子 结 
合 后 ， 虽 然 收 敛 速度 有 所 下 降 ， 但 明显 提高 了 全 局 搜索 能 力 ; 
而 本 文 算法 采用 Givens 变换 并 引入 二 阶 振荡 环节 , 即 提 高 了 收 
敛 速 度 ， 又 增强 了 全 局 搜索 能 力 ， 并 且 极 易 跳出 局 部 极 值 ， 算 
法 性 能 越 来 越 好 。 
另外 ， 相 较 与 其 他 PSO 改进 算法 ， 性 能 对 比 结果 如 图 10 
示 。 

文献 [8] 采 用 自 适 应 惯性 权重 PSO 算法 ， 全 局 搜索 能 力 较 
差 ， 且 易 陷 入 局 部 最 优 : 文献 [9] 引 入 遗传 变异 环节 后 全 局 搜索 
能 力 有 所 提高 , 但 由 于 其 复杂 度 较 高 收敛 速度 较 慢 ; 比较 之 下 ， 


Im 


Ne 
» 
| 


CT 
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本 文 算法 的 收敛 速度 有 所 提高 ,同时 具有 较 好 的 全 局 搜索 能 力 ， 


不 易 陷入 局 部 极 值 。 
30 人 
es 文献 8 改进 算法 
ol 文献 9 改进 算法 | 
一 一 本 文 改进 算法 


0 50 100 200 250 300 


迄 代 次 数 
10 ”各 算法 收敛 曲线 2 
Fig.10 Convergence curves of algorithms 
2) 分 离 性 能 
为 了 进一步 体现 本 文 算 法 的 分 离 性 能 ， 使 用 相似 系数 作为 
评价 性 能 指标 ， 即 分 离 信 号 与 源 信号 的 相似 性 。 相 似 系数 3 
为 


现 


书 


(19) 


当 》 与 9 相互 独立 时 ， 包 =0 ; 当 y》=C5; (ce 为 常数 ) 时 ， 
与 =1， 说 明 s(1) 与 y(0) 越 相似 ， 即 分 离 效 果 越 好 。 使 用 传统 
EASI 算法 、 文 献 [8] 的 自 适 应 惯性 权重 PSO 算法 、 文 献 [9] 遗 伟 
变异 算法 与 本 文 提出 的 算法 分 别 进行 仿真 ， 其 中 EASI 算法 选 
择 非 线性 函数 f(y)=y ， 将 得 到 的 平均 相似 系数 进行 比较 。 实 
验 采用 机 械 振动 信号 进行 实验 ， 结 果 如 表 1 所 示 。 

表 1 各 算法 相似 系数 对 比 


Tablel Comparison of similarity coefficients 


信号 1 平均 相似 系数 这 号 2 平均 相似 系数 

最 大 值 最 小 值 最 大 值 最 小 值 
EASI 0.9292 0.4251 0.9105 0.0973 
文献 [8] 0.9225 0.2684 0.9588 0.4007 
文献 [9] 0.8642 0.1132 0.9754 0.1433 
本 文 算法 0.9999 0.0865 0.9855 0.0009 

由 表 1 可 以 看 出 ，EASI 算法 、 文 献 [8] 算 法 、 文 献 [9] 算 法 


的 相关 系数 最 小 值 都 比较 大 ， 说 明 对 动态 信号 分 离 效 果 不 好 。 
而 本 文 提出 的 改进 算法 ， 即 基于 Givens 变换 和 二 阶 振荡 
W-C-PSO 算法 ， 其 相关 系数 的 最 大 值 为 0.9999， 相 比 于 其 他 算 
法 更 接近 1; 其 相关 系数 最 小 值 为 0.0009。 说 明 本 文 提 出 的 算 
法 具有 更 好 的 分 离 性 能 ， 为 处 理 动态 信号 提供 了 新 的 思路 和 依 
据 。 

3) 运 行 时 间 

为 了 多 方面 比较 算法 的 性 能 ， 表 2 
间 。 


全 出 了 各 算法 的 运行 时 
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表 2 算法 运行 时 间 [5] 机 志 成 , 韩 大 伟 ,， 陈雷 , 等 . 基于 复 Givens 矩阵 与 蝙蝠 优化 的 卷 积 讶 分 
Table 2 Algorithm run time 离 算法 [四 , 通信 学 报 , 2016, 37 (7): 107-117. (Jia Zhicheng, Han Dawei, 
算法 运行 时 间 /s Chen Lei, et al. Convolutive blind separation algorithm based on complex 
W-C-PSO 12.565630 Givens matrix and bat optimization [J]. Journal on Communications, 2016, 
二 阶 振荡 +W-C-PSO 12.703810 37 (7): 107-117. ) 
本 文 算法 12.385740 [6] 陈雷 ,， 甘 士 忠 ， 张 立 毅 ,等 . 基于 样 条 插值 与 人 工蜂 群 优化 的 非 线性 育 
表 2 看 出 ， 各 个 算法 的 运行 时 间 很 接近 ， 虽 然 引入 二 阶 源 分 离 算法 [J]. 通信 学 报 ，2017，38 (7): 36-46. (Chen Lei，Gan 


振荡 粒子 使 得 算法 的 收敛 速度 较 W-C-PSO 算法 变 慢 ， 但 是 由 
图 9 可 知 W-C-PSO 算法 的 性 能 并 不 高 ， 所 以 本 文 算法 在 此 基 
础 上 使 用 Givens 矩阵 表示 分 离 矩阵 ， 由 3.1 节 可 知 ， 本 文 算法 
减少 了 n(n+])/2 个 参数 ， 使 得 运行 时 间 有 所 减少 ， 加 快 了 收 


5 ”结束 语 


本 文 在 传统 PSO 算法 上 对 惯性 权重 和 学 习 因 子 进行 结合 ， 


增强 算法 的 整体 性 ， 提 高 了 算法 的 全 局 搜索 能 力 ;， 在 此 基础 上 
又 引入 二 阶 振荡 环节 ， 增 加 了 粒子 的 多 样 性 ， 使 其 不 易 陷入 局 


部 极 值 ; 并 采用 Givens 变换 降低 了 算法 的 复杂 度 , 提高 了 算法 
收敛 速度 。 仿 真 表明 ， 对 机 械 振动 信号 和 语音 信号 盲 分 离 均 具 
有 较 好 的 效果 ; 在 收敛 性 、 全 局 搜索 能 力 以 及 分 离 性 能 等 方 
本 文 算法 表现 更 优 ， se 
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